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The Supernova Legacy Survey (SNLS) has produced a high-quality, homoge- 
neous sample of Type la supernovae (SNe la) out to redshifts greater than z = 1. 
In its first four years of full operation (to June 2007), the SNLS discovered more 
than 3000 transient candidates, 373 of which have been confirmed spectroscopi- 
cally as SNe la. Use of these SNe la in precision cosmology critically depends on 
an analysis of the observational biases incurred in the SNLS survey due to the in- 
complete sampling of the underlying SN la population. This paper describes our 
real-time supernova detection and analysis procedures, and uses detailed Monte 
Carlo simulations to examine the effects of Malmquist bias and spectroscopic 
sampling. Such sampling effects are found to become apparent at 2; ~ 0.6, with 
a significant shift in the average magnitude of the spectroscopically confirmed 
SN la sample towards brighter values for z > 0.75. We describe our approach 
to correct for these selection biases in our three- year SNLS cosmological analysis 
(SNLS3), and present a breakdown of the systematic uncertainties involved. 

Subject headings: cosmology: observations — supernovae: general — methods: 
data analysis — techniques: photometric — surveys 



Introduction 



Type la supernovae (SNe la) have frequently demonstrated their value as standardizable 
candles out to redshifts of 2; > 1.5. Stu dies of the expansion rate of the universe u sing SNe la 
have revealed the cosmic acceleration (jPerlmutter et al.l Il999l : iRiess et al.l Il998l ) , providing 
the most direct evidence for the presence of Dark Energy needed to drive the expansion of 
the u niverse. Subsequent searches have led to SN la discoveries over a wide range of redshifts 



le.g- 



Astier et al.ll2006l : Miknaitis et al.l 120071 : ISako et al.ll2008l : iRiess et al.l 120071 ). 



SNe la are relatively rare, transient objects lasting a few weeks in the rest-frame. They 
are most efficiently found through photometric monitoring of large areas of sky using a 
"rolling- search" tec hnique to both find and follow SN events. Cur rent examples include the 
ESSENCE survey jMiknaitis et al.ll2007l : IWood-Vasev et al.ll2007h . the SDSS-I I SN Survey 



(|Sako et al.ll2008l : iKessler et al.ll2009l ). and the Supernova Legacy Survey (SNLS: lAstier et al. 
20061 ) . The SNLS observed four 1- square- degree Deep Fields as part of Canada- France- 



Hawaii Telescope Legacy Survey (CFHT-LS)0. Multi-color light curves were obtained using 
the MegaCam imager on CFHT, with follow-up spectroscopy from the Gemini, VLT, and 



^http://www. cfht.hawaii.edu/Science/CFHLS/ 
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Keck telescopes. Owing in large part to its careful survey design and a substantial com- 
mitment of observing time from CFHT and the follow-up telescopes, the SNLS now offers a 
large, consistent, and well-defined set of SNe la measured out to 2; ~ 1. 

Real-time supernova searches require adequate light-curve sampling, a rapid output of 



up studies (e.g., ISuUivan et al.l l2006al : iMiknaitis et al. 



advantageous to obtain spectroscopy near maximum SN brightness (e. 



. 2008 


. It is most 


g., see 


Howe 


et al. 


Conlev et al. 


2006 





20051 ). Since SNe la have typical rise times of ~ 19 days (rest-frame B, IConley et al.l 
and considering that observations are typically only taken during dark and grey time, in 
some cases only a couple of measurable epochs are available before peak. Data organization 
and the coordination of follow-up activities are therefore critical. In the SNLS, Canadian and 
French teams each ran independent detection pipelines, then amalgamated both candidate 
lists into a centralized database. Herein, we report exclusively on the Canadian real-time 
analysis pipeline. 

As the extent and quality of SN la data sets improve, the relative impact of statistical 
uncertainties decrea s es and the management of systematic error s becomes the limiting factor 



(IKessler et al.l l2009l : iRegnault et al.l l2009l : IConley et al.l I2OIOI ). One important source of 



systematic error in SN la studies arises from the incomplete sampling of the underlying 
SN la population: surveys tend to preferen tially discover an d follow brighter objects at 
higher redshifts, leading to Malmquist bias (iMalmquistI Il936l ) and other selection effects. 
Using detailed Monte Carlo simulations of artificial SNe la in all of the real-time detection 
images, we can determin e the recovery effic iency and calculate the magnitude of the various 
sampling biases (see also lConley et al.ll2010l ). These can then either be directly applied to the 
data, or used to motivate fitting pri ors or corrections used when determining cosmological 
parameters (e.g.. iKessler et al.ll2009l ). 



This paper is organized as follows: An overview of the SNLS and a description of the 
Canadian SNLS real-time-analysis (RTA) pipeline and web database are given in ^ ^ 
presents the methods used in ranking candidates for spectroscopic follow-up observations. 
The Monte Carlo simulations used to calculate the SNLS detection efficiencies are described 
in §1] along with the recovery results. The effects of Malmquist bias and spectroscopic 
sampling are explored in ^ with a discussion of the systematic errors in the Malmquist bias 
in ^ The paper concludes with a summary in ^ 
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2. Observations and real-time analysis 



2.1. SNLS observations 



The five-year Supernova Legacy Survey ran from 2003-2008 and was based on images 
obtained from tlie CFHT-LS Deep Synoptic Survey. Tlie four Deep Field pointings were 
positioned to avoid extremely bright stars, minimize Milky Way extinction, and overlap with 
the fields of other multi-wavelength surveys to provide complementary data. The coordinates 
of the SNLS Deep Field centers are given in Table [TJ 

Observations were obtained in qu eue-scheduled mode at CFHT using the MegaCam 
wide-field imager feoulade et aPboOsh . MegaCam consists of 36 2048 x 4612 pixeP CCDs 
arranged in a 4 x 9 square mosaic. The CCDs have a resolution of ~ 0.186"/pix and provide 
good sampling of the 0.7" median seeing at CFHT. A large (1.5') dithering pattern was used 
to cover the 80" wide gaps between CCD rows in the MegaCam mosaic. 

The SNLS employed a rolling-search observing strategy, with target fields imaged every 
3 — 4 days in dark or grey time using multiple filters to provide optimal light-curve sampling 
and the opportunity for early SN discovery. This cadence translates to ~ 2 — 3 days in the 
SN rest frame. Each MegaCam queue run lasted on average 16 — 17 nights centered around 
the new moon, typically providing 5 epochs (nights) of observation per field per lunation. 
An advantage of this scheme is that the same telescope is used to discover and follow-up all 
candidates, leading to considerable advantages in calibration. 

The four filters used by SNLS — gurMHiZM — are similar (but not identical) to the 



filter system adopted by the Sloan Digital Sky Survey flFukugita et al.l Il996l : ISmith et al. 



2OO2I ). Deep Field images were also obtained in um, but these were not time-sequenced and 
thus were not used in the real-time analys i s. A d etailed description of the SNLS photometric 
calibration is provided in iRegnault et al.l ( 120091 ). Observations in the primary search filter 



Table 1. Deep Field centers 



Field 


RA (J2000) 


DEC (J2000) 


Overlapping Surveys 


Dl 


02:26:00.00 


-04:30:00.0 


VIMOS, SWIRE, GALEX, XMM Deep 


D2 


10:00:28.60 


02:12:21.0 


COSMOS/ACS, VIMOS, SIRTF, GALEX, XMM 


D3 


14:19:28.01 


52:40:41.0 


(Groth Strip) ACS, DEEP-II, GALEX, AEGIS 


D4 


22:15:31.67 


-17:44:05.7 


XMM Deep 
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[hi) were alternated between "long" and "short" epochs in order to balance the requirements 
for adequate light-curve coverage with the need to conduct several deep searches within each 
queue run. Filter spacing and nightly exposure times for a typical queue run are provided 
in Table H 

A major strength of the SNLS program is that it obtained better filter coverage and light- 
curve cadence than had been achieved by previous high- 2; surveys. Other advantages included 
having a large quantity of dedicated telescope time with data acquisition coordinated through 
CFHT's Queued Service Observing mode. Furthermore, developing a careful strategy for the 
organization, analysis, and communication of SNLS data products was necessary to efficiently 
manage the continual acquisition of new target information. The process employed by the 
SNLS collaboration included linking remote search pipelines with a centralized database, 
and providing web-based applications as a simple user interface with communication and 
analysis functionalities. These are described in the next section. 



2.2. Image processing 



MegaCam images of the Deep Fields destined for SNLS analysis w ere pre-processed in 
real-time using the CFHT Elixir pipeline ( IMagnier fc Cuillandrd I2004J ) . The Elixir pipeline 
performed bias subtraction and flat-fielding, as well as a basic fringe subtraction on the 
and zm images. Elixir-processed images were promptly made available to the SNLS team 
for further processing. 

The real-time detection pipeline implemented by the Canadian SNLS team is outlined in 
the flowchart of Figured! Many of the processing steps were accomplished using IRAF@ tasks. 
Due to time constraints and observatory requirements, the image processing was carried out 
remotely on SNLS computers at the CFHT headquarters in Waimea. Certain shortcuts were 
taken in the real-time pipeline as comp ared with the "fin al" photometry procedure used in 
the cosmological analysis (described in iGuy et al.l |2010| ). These were necessary since the 
SNLS RTA was operating under significant time pressures for discoveries and analysis. 

Elixir-processed images were trimmed and reoriented, then bad pixel masks were created 
to flag and remove image artifacts and saturated stars. To account for optical distortions, 
astrometric solutions were fit to each CCD using combined catalogs from USNO-B and SDSS 



^Image Reduction and Analysis Facility (IRAF) is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under 
cooperative agreement with the National Science Foundation. 
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MegaCam Raw 
Images (36 CCDs) 




CFHT Elixir 
Processing 

Flatfielding, 
fringe correction 




5-1 images aligned to 
ICRS/SDSS to +/-0.06" 



SWARP 

Geometric resampling, 
coadd images, 
CR rejection 



I 



Elixir-processed 
Images (36 CCDs) 



Image 
Preparation 

Trim Si reorient CCDs, 
mal<e bad-pixel masks 



Astrometric 
Calibration 

Fit WCS using 
internal star catalogs 



Pixel Area 
Correction 

Flux conservation 
(distortion) 



Flux Scaling 

Secondary standards 
(to common epoch) 



PSF Match 

Match PSFs of 
epoch and 
reference images 



Difference 

Subtract current 
epoch - ref = diff 



I 



Photometry 
Pipeline 

Measure PSF and 
aperture fluxes 



Detfinder 

Search for objects 
in diff, cull out junk 



Photometry 
for all 
objects 
in current 
epoch/filter 



Detection 
Review 

Create mosaics of 
detection "triplets" 



Visual 
confirmation 
of candidates 



EPOCH 




List of new 
candidates 
in current 
epoch 



REF 



DIFF 



Fig. 1. — The Canadian RTA pipeline. The dotted box dehneates the procedures that 
monitor incoming Ehxir-processed images. The rest of the steps proceed on the pre-processed 
images as described in the text. 
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to align the images to the International Celestial Reference System (ICRS). The resulting 
fits had RMS residuals that were typically better than < 0.06" (0.3 pixels) across the mosaic. 

Photometric alignment of the RTA images was accomplished using aperture flux mea- 
surements of ~ 20 000 — 40 000 "RTA secondary standards" selected over a wide range in 
brightness for all filters within each field. These secondaries were chosen from isolated objects 
exhibiting no significant changes in brightness over many epochs. Comparing the measured 
fluxes for the secondaries in a given epoch to their catalog values provided a multiphcative 
scaling factor used to flux calibrate each image. As such, all flux measurements over the 
duration of any SN light curve were consistently in the same reference system, removing vari- 
ations in exposure time, airma ss, and extinction. In the RTA, the images were calibrated 



to the AB photometric system ( Oke fc Gunnll983); th is differs from the me thod used in fi- 



nal p hotometry, which calibrated to the iLandoltl ( 1l992l ) photometric system (IRegnault et al. 



20091). 



Individual fiux-scaled CCD images were adjusted for variations in effective pixel area 
using the IRAF MSCRED (mosaic data reduction) task MSCPIXAREA. Resampling and co- 
addition of the night's images were then done using the TERAPIX SWarp|f| program to form 
a median mosaic image and its accompanying weight-map for each field-|-filter combination. 

The subtraction of a deep reference image from the science mosaic revealed the presence 
of SNe and other variable objects in the data. References for each field were constructed 
by selecting photometric images with excellent seeing (< 0.5") and good transparency. To 
avoid the presence of SN light in the reference images, only frames taken in years prior to 
the current observations were included in each stacl<0. These component images were flux 
scaled, WCS-aligned and co-added as described above to form the reference image. During 
the course of the survey, new RTA stacks were constructed at least once per season as more 
images became available to improve the depth of the reference mosaic. 

Differences in image quality between the science and reference mosaics must be taken 
into account prior to image subtraction. A non-parametric point-spread function (PSF) 
matching routine was used to effectively degrade the seeing of the bette r quality image 



typically the reference — to match the seeing of the other image (see iFraser et al.l 12008 
Appendix A). A Moffat function was fit to a select set of a few thousand isolated, bright 
"PSF stars" , culling those with poor-quality fits. A variable convolution kernel was computed 



" http : //astromat ic . net 

^RTA detections made during the first year targeting each field were made using whatever prior data was 
available. The quality of the reference stacks improved with each successive observing season. 
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for each star, matching the two PSFs. The reference image was then convolved with the 
kernel to produce a PSF-matched mosaic. Finally, the matched frame was subtracted from 
the science frame. 

A sample "detection triplet" is presented in Figure [2] for a high-redshift SN candidate: 
SNLS-04D2ca at z = 0.835. It shows thumbnail images of corresponding regions in the 
science image for one epoch, the PSF-matched reference image, and the difference. 

One drawback of the real-time analysis procedure is that, in the interest of processing 
speed, the method of co-adding the individual epoch images prior to the PSF-matching step 
meant that if the image quality changed significantly during a night's set of observations, 
the resulting difference image may have contained residuals from poor subtractions of host 
galaxies. Normally, such residuals were easy to identify and immediately cull from the output 
detection list, as described in the next section. 



2.3. Supernova detection 

Once the host galaxy light was removed via the subtraction of the reference image, 
potential supernova candidates appeared in the subtracted images as star-like objects with 
positive flux (see the right-hand panel of Figure [2]). The subtracted im ages were searched 



using an algorithm that employed the TERAPIX SExtractor routine (IBertin fc Arnouts 



19961 ) to build a catalog of possible objects and apply a series of tunable culling procedures. 
Each preliminary catalog of detections was cleaned of objects with very low signal-to-noise or 
a significant fraction of bad pixels, as well as objects with irregular profiles or high second- 
order moments. Adjustable image quality limits were available for additional culling on 
PSF-fit quality parameters and background sky measurements. 

Residual features in the difference images could result from poorly subtracted galaxies 
or bright stars. These residuals commonly appeared as alternating bright and dark regions 
— the bright areas could sometimes survive the other culling procedures. To compensate for 
this, the detection code also compared the image statistics in defined regions immediately 
surrounding the candidate to help remove such false detections from the object list. 

The resulting raw list of variable object candidates in each epoch was then used to 
produce a composite mosaic of detection triplets for visual review and a manual selection 
of candidates. The number of raw output detections depended on the parameter limits 
established for the cuts (e.g., the S/N cutoff). We generally opted to use liberal cuts and 
do a manual pre-selection of SN candidates from the resulting < 200 raw detections. Of 
these, typically ~ 10% remained as plausible SN candidates. Using the visual display of 
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Table 2. Typical SNLS/Deep exposures 



Filter 


ExpTime ( 


sec) 


Epoch Number 


Qm 


5 X 225 = 


1125 


1,2,3,4,5- 




5 X 300 = 


1500 


1,2,3,4,5 


iM (long) 


7 X 520 = 


3640 


1,3,5 


iu (short) 


5 X 360 = 


1800 


2,4 




10 X 360 = 


3600 


1,3,5 



^qm acquisition depended on Moon brightness. 




Fig. 2. — A sample detection triplet showing an SN la candidate (SNLS-04D2ca) aX z = 
0.835. The science image in the im filter for the night of 2004-03-10 is shown on the left, 
the PSF-matched im reference image is in the center, and the difference between the two is 
the image on the right. The subtracted mosaic is used for candidate detection as described 
in the text. Each image is 19" wide. 
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candidates, many of the remaining PSF-matching residuals, image artifacts, moving objects, 
and stellar cores that may have survived the previous culling steps could be removed. 

Detections were primarily made using the im (and to a lesser extent tm) images in 
order to focus on the target redshift range of interest. Candidates that were too faint in tM 
were poorly-suited for SNLS follow-up spectroscopy and cosmological analysis. A detailed 
analysis of the detection incompleteness is presented later in ^ 

2.4. Photometry 

Supernova flux measurements were calculated using a PSF-fitting routine. The same 
sets of isolated stars used in PSF-matching step ( §2.2p were used to create a mean local 
PSF and residuals map for input to the fitting routine. The stellar model was then scaled 
to match each detection in the difference image, thus providing the flux of each source. In 
addition to the PSF photometry, aperture fluxes were also calculated using a range of static 
and seeing- dependent aperture sizes. 

The output PSF-fit fluxes measured in each epoch were used to construct multi-color 
light curves for all of the variable-object candidates detected by SNLS. These light curves 
were then fit with SN la templates to assist in classification and to provide photometric 
redshift measurements; this process is described in ^ 

A sample set of SN la light curves in gui^MiM is shown in Figure [3l demonstrating the 
typical observing cadence obtained by the SNLS. 

2.5. The SNLS web database 

Once the manual review and culling was complete for a given epoch, the preliminary 
detection list was loaded into the SNLS database in preparation for ranking and follow-up 
measurements. New SN candidates and photometry were made available to the community 
via the World Wide Wet|f|, typically within 4 — 6 hours per field. 

An observing request utility was used communicate the latest SNLS requirements to the 
CFHT queued service observers (QSOs). This tool kept the QSOs advised of the detailed 
priorities of the survey, reaching a level beyond what could normally be expressed in the 
Phase II application that was submitted at the start of each run. Often, observing condi- 



" http : //legacy . astro . utoronto . caj 
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Date 



Sep Oct Nov Dec Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec Jan 




Sep Oct Nov Dec Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec Jan 



2005 



Fig. 3. — Observed light curves for some of the SNLS SN la candidates, showing tM (top), 
tm (middle), and qm (bottom). Each observation is given by a dot on the corresponding 
SN la light-curve fits (solid lines). The observations were clustered during dark time, with 
additional gaps due to runs of poor weather. 



- 12 - 



tions, SN candidate status, and follow-up telescope availability called for some flexibility in 
field/filter observing priorities for SNLS. 

Additional online tools facilitated the coordination of observing schedules and follow- 
up allocations so as to optimize the use of the available 8-10 m telescope time. A detec- 
tion "forum" was used to collect comments and personalized rankings for candidates being 
considered for follow-up spectroscopy. Automated (but customizable) online finder charts 
simplified subsequent observations in both single-target and multi-object modes. 

It is worth noting that the light curves of all candidates in the database were followed 
until they became too faint or were unequivocally rejected as a non-SNe. No candidates were 
thrown away, no matter how poor; all have been stored in the database for later studies. 



3. Prioritization and spectroscopy 

As with most supernova surveys, the SNLS was limited by the amount of time available 
for spectroscopic follow-up. Even with its maximum allocation of ~ 140 hours of 8-10 m 
telescope time per semester, the SNLS was only able to follow a small fraction of all possible 
SN detections. Differentiating the SNe la from other variable sources using only a few early 
observations presented a significant challenge. 



3.1. Candidate ranking 



The success of any supernova survey relies on the careful selection of SN la candidates 
with enough time to optimally schedule spectroscopic follow-up close to peak magnitude. 
Light-curve fits and brightness cuts were the primary methods used to identify and rank 
SNLS targets for spectroscopy. The addition of photometric pre-selection dramatically im- 
proved the SNLS s pectroscopic efficiency over comparable surveys, even at high redshift 
jHowell et allboosh . 



Probable SNe la could typically be identified from two or more epochs of real-time pho 



tometry with the supernova light-curve fitting technique described in detail in ISuUivan et al. 



(l2006al ). The output light-curve fits are parameterized by redshift, stretch, peak brightness, 
and time of maximum light. Compar isons between photorn etric and spectroscopic re dshifts 
for the SNLS SNe la are available in ISullivan et al.l ( l2006al ) and lPerrett et al.l (120101 ). Non- 
SNe la could be screened out based on their measured colors {qm — '^m) and rise times. 
All SN candidates were measured against several criteria to determine their suitability and 
prioritization for spectroscopic follow-up. Candidates with poor-quality SN la light-curve 
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fits, or with fitted stretch values outside the range 0.7 < s < 1.3, were generally given the 
lowest ranking or were rejected entirely. 

Observing near maximum light yields the brightest contrast between the SN and its 
host, and is the time when the spectral lines most clearly differ between supernova types. 
SNLS SNe la were normally discovered before peak, and the predicted phase from the fits 

allowed for optimal spectroscopic scheduling. Each candidate SN was required to have been 
detected at most seven days (rest frame) past its estimated peak date in order to qualify for 
timely spectroscopic observing. Objects discovered at the very start or end of an observing 
season were typically discarded if they would have no measurements either before or after 
peak. Enough of the candidate's light curve must be measured to provide adequate fits and 
prove useful for cosmology. 

Limits were placed on the minimum observed brightness and fractional increase of each 
candidate compared with its host galaxy. Percent increase is defined as: 

%inc = l^SL^Jl^ X 100%, (1) 

/ref 

where flux is measured in small (3-pixel radius) apertures in the science image (/sci) and the 
reference (/ref)- 

Due to the vagaries of weather and telescope scheduhng, it is impossible to give a 
precise description of the selection process for sending candidates for spectroscopic follow- 
up. However, a good description of the average algorithm is as follows, with all observed 
magnitudes presented in the AB system: 

1. Extremely bright SNe (zm < 22.9) always qualified for spectroscopic follow-up. 

2. Moderately bright SNe (22.9 < in < 23.8) quahfied if they had %inc> 30 in some 
epoch. 

3. Fainter SNe (23.8 < iM < 24.4) must have had %inc> 100 in some epoch to quahfy. 

Candidates fainter than Im ~ 24.4 were generally never observed spectroscopically. 

The basic form of these hmits was suggested by the experience of those selecting can- 
didates for follow-up, and tuned by comparing the distributions of peak magnitude, color, 
and light-curve width of moderate redshift SNLS SNe {z < 0.6), where the survey should be 
essentially complete, to those of higher redshift SNe. 

In practice, the selection process remained somewhat subjective. External factors also 
played a role in the ranking of candidates by SNLS members — e.g., the availability of 
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follow-up telescope time, detailed scheduling limitations, long-term weather conditions, and 
the desired spread in redshifts needed for cosmological measurements. 



3.2. Spectroscopy 

Spectroscopic observations are necessary to verify candidate types and to obtain accurate 
redshift measurements for the SN and/or i ts host galaxy. SNLS spectra were obta ined from 



8-10 m class telescopes in cluding Geminj§ (Howell et al. 



2005 



iBronder et a. 



Large Telescope (VLTjil fealland et aPboOoh . and Kec#l JeUIs et al.lboosl h 



20081), the Very 



The redshift distribution of the spectroscopically-identified SNe la found by SNLS is 
shown in Figure |H The spectroscopic sample is fairly complete out to 2; ~ 0.6, at which 
point the numbers begin to drop off as the amount of exposure time required per candidate 
increases significantly. This effect on the incompleteness of the SNLS is examined later in 



Histograms showing the breakdown of the sample by type and spectroscopic redshift 
from the various telescopes are presented in Figure [51 Candidates labeled as "SN?" are 
those that have not been confirmed as SNe la or core-collapse SNe; this category may also 
include other types of variable objects that were not easily distinguished from supernovae. 
The non-SNe category mostly includes active galactic nuclei (AGNs) and other objects, most 
of which were observed prior to the implementation of the photometric fitting technique for 
follow-up ranking described in §3.11 



The median spectroscopic redshifts of the SNe la observed at the three observatories are 



^This work is based in part on observations obtained at the Gemini Observatory, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on 
behalf of the Gemini partnership: the National Science Foundation (United States), the Science and Tech- 
nology Facilities Council (United Kingdom), the National Research Council (Canada), CONICYT (Chile), 
the Australian Research Council (Australia), CNPq (Brazil) and CONICET (Argentina). Gemini pro- 
gram IDs: GS-2003B-Q-8, GN-2003B-Q-9, GS-2004A-Q-11, GN-2004A-Q-19, GS-2004B-Q-31, GN-2004B- 
Q-16, GS-2005A-Q-11, GN-2005A-Q-11, GS-2005B-Q-6, GN-2005B-Q-7, GN-2006A-Q-7, GN-2006B-Q-10, 
and GN-2007A-Q-8. 

^Observations made with ESO Telescopes at the Paranal Observatory under program IDs 171.A-0486 
and 176.A-0589. 

^Some of the data presented herein were obtained at the W.M. Keck Observatory, which is operated 
as a scientific partnership among the California Institute of Technology, the University of California and 
the National Aeronautics and Space Administration. The Observatory was made possible by the generous 
financial support of the W.M. Keck Foundation. 
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Redshift 

Fig. 4. — The number of spectroscopic SNe la from the SNLS as a function of spectroscopic 
redshift. This plot includes confirmed SNe la (shaded region) and probable SNe la (open 
region) discovered up to the end of D3 observations in June 2007. The solid curve represents 
the rise in number expected purely due to the increasing volume at greater redshifts, scaled 
to match the observed value in the z = 0.5 — 0.6 bin. 
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Fig. 5. — Supernova candidates from the SNLS with spectroscopic follow-up observations. 
Top: Candidate types as determined by spectroscopy from Gemini, VLT and Keck. Bottom: 
Redshift distribution for Type la SNe from the SNLS sample, as in Figure H] but separated 
by telescope. Both plots only include discoveries up to the end of D3 observations in June 
2007. The SN la category includes the probable (SN la?) candidates. 
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z = 0.72 (Gemini), z = 0.62 (VLT), and z = 0.52 (Keck). Higher-redshift candidates were 
more frequently sent to Gemini because of the larger observing overheads, and because the 
nod- and- shuffle mode available at Gemini is useful for removing sky lines to obtain better 
spectra of high- 2; SNe. The Keck observations included several SN la programs with different 
science goals. 



3.3. Candidate numbers 



In its first four years of full operation (to July, 2007)|^, the SNLS discovered > 3000 
astrophysical transient events. Of these, 373 were spectroscopic ally identified as SN e la. 



This is the same SNLS sample used in the SN la rates analysis of iPerrett et al.l ( I2OIOI ). As 
expected, there is a strong correlation between the number of SN la detections and the total 
observing time in iM, as demonstrated in Figure [61 This figure contains four plots, one per 
field, each showing the running total of SNLS exposure time per filter (with {m in red) for 
comparison with the cumulative number of SNe la discovered (the shaded areas). 

A summary by type for all SNLS transient detections (excluding moving objects and 
subtraction artifacts) is provided in Table |3l A^ refers to the number of objects of each type in 
the SNLS real-time database and A^spec is the number of those with available spectroscopy. 
Trailing question marks "(?)" in the type label in Column 1 indicate that the number 



^In June 2007, the iM filter for MegaCam was destroyed during a malfunction of the filter jukebox. 
Candidates discovered after this date were observed with a new im filter, requiring new references and 
calibrations, and are still currently under analysis. Therefore, they are not included in the sample described 
in this paper. 



Table 3. SNLS candidates to 2007-06-30 



Type 


N 


A^ 

^ ' spec 


SN Ia(?) 


373 


373 


SN II(?) 


51 


50 


SN Ib/c(?) 


17 


17 


Other SN(?) 


1309 


166 


non-SN 


1360 


22 


Total 


3110 


628 



- 18 - 




Fig. 6. — A timeline showing the cumulative number of probable and confirmed SNe la 
discovered in each Deep Field by SNLS (shaded area), along with the total exposure time 
obtained in the separate filters. 
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includes both confirmed and probable identifications of that type. SNe la? are classified 
as "probable" because their spectra best match that of a Type la, but another type of 
supernova — typically a SN Ic — cannot be conclusi vely ruled out from the s pectrum alone. 
SNe l a? are equivalent to the SNe la* described in iHowell et al.l (|2005[ ) and iBalland et al. 
( I2OO9I ). New SN-like detections were given the default label "SN?" in the database until 
they could be identified via spectroscopy. (The one exception is a SN II-P SN, which could 
be identified from the obvious plateau in its light curve.) Three of the objects categorized as 
SN Ib/c(?) are actua lly confirmed SNe Ic: SNLS-03D4aa, SNLS-04D4jv, and SNLS-06D4dv 
tealland et aPbopgh . 



We believe the bulk of the "Other SN" category to be real SN events that were not 
followed spectroscopically due to a lack of follow-up time. Many of these are very faint (in 
some cases below the spectroscopic limit) and likely high redshift; some had colors that were 
clearly incompatible with SNe la. Others will be SN la events in our redshift range. These 



candidates will be analyzed further to measure the volumetric SN la rate in iPerrett et al. 



( I2OIOI ). The "non-SN" category includes a large number of other astrophysical transient 
events, including AGN and variable stars, identified via their erratic light curves and central 
location in their "hosts" . 

We emphasize that the relative numbers of the different types of SNe in Table [3] are 
strongly affected by spectroscopic selection procedures. The total numbers are also influenced 
by factors such as losses due to poor winter weather conditions on Mauna Kea, as well as 
the slow ramp-up of observing during the pre-survey period. These observational effects are 
incorporated into the candidate recovery statistics that are measured in the next section. 



4. SNLS Detection Efficiency 

A thorough understanding of the recovery statistics in the RTA pipeline provides the ba- 
sis for calculating sample selection biases. Detailed Monte Carlo simulations were carried out 
to measure the efficiency of SN la recovery from the SNLS real-time images. Although these 
simul ations could have b een performed concurrently with the real-time candidate searching 



(e.g. 



Dilday et al.l 120081 ). in the SNLS we perform these simulations completely separately 
from the daily search activities. A description of the input test population is provided in 
§4.H and the results of the simulations are outlined afterwards in §4.21 
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4.1. SNe la input parameters 



In all, 15 realizations of 10 000 artificial SNe each were performed individually for every 
Deep Field observing season up to the end of D3 in the summer of 2007. A "field-season" 
represents the period over a year during which a given field was visible: in this current 
sample, D3 was visible for five seasons (counting the pre-survey observing), while the rest of 
the fields were each visible for four. This yields a total of more than 2.5 x 10^ artificial SNe 
tested in the simulations during 17 field-seasons. 

To begin, a random redshift was selected for each input artificial SN from a uniform 
distribution over 0.1 < z < 1.2. The artificial SNe were then randomly assigned to hosts 
chos en from the SNLS fiel d galaxy catalogs compiled from measurements of deep image stacks 
Sullivan et al.ll2006bl ). The limiting magnitude of the reference stacks is ijv/ ~ 25.7 mag 



[see 

(AB), much fainter than the detection limit for any given epoch. Photometric redshifts for 
the field galaxies were calculated based on the PEGASE.2 galaxy spectral evolution code 
( iLe Borgne et al.l 120041 ). A random host galaxy within 2; ± 0.02 of the selected redshift 
was picked in the appropriate Deep Field, weighted by a probability approximated by the 
two-component model for SN la rates (SNR): 



P{z) ~ SNR{z) = Ax M{z) + Bx SFR{z) 



(2) 



f lMannucci et al.l 120051 : IScannapieco fc BildstenI 120051 ). A = 5.3 x 10"^^ SNeyr"^ Mg^ is the 
mass coefficient, and B = 3.9 x lO'^SNe yr"^ (M^ yr ^) ^ is the coefficient of the star- 
formation rate (SFR) term (adopted from ISullivan et al.ll2006al ). As a precaution against 
overcrowding, a field galaxy was not permitted to host more than one artificial SN during 
any given season. 

The artificial SNe were assigned galactocentric positions drawn from a 2D Gaussian 
distribution about the host centroid^. verifying that each fell within the galaxy's detection 
ellipse returned by SExtractor. This produced a list of input coordinates for the input sample 
according to the astrometric solution of the reference stack. 

The input objects were then randomly allocated stretch values drawn from a uniform 
distribution in the range 0.5 < s < 1.3. This wide range was adopted in order to over-sample 
the parameter space of interest, to include recovery statistics for potential sub-luminous SNe 



^°The galaxy light in the stacks has already been convolved with the seeing (~ 0.5"), yet the SN positions 
are accurate to ^ 0.1". This will slightly underestimate the true central concentration of the input SN 
distribution. However, we find that the recovery fr action is quite flat w ith radius, as shown later in Figure mi 
This is also in agreement with the conclusions of Howell et al. ( 2000l ) , who demonstrate that selection bias 
at small host radii is not significant for CCD observations of high-z SNe la. 
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in the SNLS data (iGonzalez-Gaitan et al.ll2010[ ). This population could later be resampled 



to investigate the observed distribution ( §4.2p . 

Each candidate was assigned a random peak date within the range min(MJDobs) — 20 < 
MJDpeak < max(MJDobs) + 10, where MJDobs are the observed dates of all of the SNLS 
images obtained for the field and year being tested. The I20 values are included to account 
for candidates that were detected while fading at the start or rising at the end of a field's 
observing period. Observational culls based on a required minimum number of observations 
before and after peak could be applied later as needed to test for selection bias. 

The peak B magnitude appropriate for a Type la at each chosen redshift was calculated 
by incorporating the stretch-luminosity relationship and an absolute magnitude derived from 
cosmology fits assuming Hq = 70. The adopted cosmological parameters are Qm = 0.27, 
Qx = 0.73, and = 0. A random dispersion in the peak 5-band magnitude (Am^ = Amag) 
was then applied after the stretch correction, drawing from a Gaussian distribution with an 
intrinsic scatter of Uint = 0.15. Again, a broad dispersion was adopted so that the output 
distribution could be resampled as needed. 



A rest-frame spectrum was extracted from an SN la spectral template (IHsiao et al. 



20071 ) and was scaled up or down to the appropriate B magnitude for each artificial object 
(with dispersion applied), thereby also affecting U and V. The output peak B magnitude 
was found by integrating over the Landolt filter functions. The U — B and B — V colors at 
peak we re calculated using a b ilinear model of the stretch-color relationship with a break at 



0.8 (iGarnavich et all 12004 ) . such that: 



-0.039 (s > 0.8) 

B-V= i V - ; 3 

[1.711 -2.187s (s<0.8) 
This s — c relationship model was derived from fits of color as a fun ction of stretch for the 



sample of low-z SNe la and the SNLS spectroscopic sample at z < 0.6 (iGonzalez-Gaitan et al. 



2OIOI). 



This procedure yields U and V rest magnitudes, and thus raw colors, for the artificial 
SNe. Gaussian noise with a dispersion of cr = 0.04 was added independently to U and V 
to a djust the colors, a nd the spectrum was "mangled" (i.e., color-corrected) to those values 



see 



Hsiao et all 120071 ) 



The spectral energy distribution was next reddened to account for Milky Way extinction, 
E{B — V)mw with -Rv,MW = 3.1, by randomly selecting from a one-sided (always positive) 
Gaussian distribution with a = 0.015. Host extinction, E{B — V)host, drawn from a uniform 
color distribution within the range —0.25 < c < 0.5, was also incorporated by assuming 
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-Rv,host = P — 1-0, with j3 = 3.0 to match the SNLS spectroscopic sample at low redshift. 
Note that th is does not preci s ely match our curr ent understanding of the SN la color-color 



relationship ( IGuy et al.l 120071 : IConley et al.l 120081 ) . which does not entirely match any mea- 



sured dust law. However, this has a negligible effect on our estimates of Malmquist bias or 
SN rates as discussed later. 

Finally, realistic SN la light curves were generated using a similar process as described 
above for calculating the peak magnitudes and colors. Randomly selected, isolated PSF 
stars in each object's CCD (brighter than the object to add) were used as the artificial 
SN models. Poisson noise was added to account for image gain. The stellar models were 
scaled and inserted with their appropriate phases, magnitudes, and positions into every 
pre-processed im image obtained by SNLS. Due to the transient nature of SNe la, most of 
the input objects added to any particular image were fainter than the detection limit, and 
no overcrowding effects were observed in the trials. The images were then processed and 
analyzed using the same RTA procedure as described in §2.21 The output detection lists for 
each epoch were compared with the input records to calculate the recovery efficiencies as a 
function of several supernova parameters. 

The strength of this straightforward yet rigorous approach is that it makes use of all of 
the actual images obtainedhj SNLS to measure recovery rates. This avoids the need to make 
assumptions about the detailed frequency and spacing of the observations, technical losses 
(e.g., missing CCDs), weather conditions, other observational effects, or to use simplified 
parameterizations to represent the recovery efficiencies. 

Before calculating the recovery statistics, the input distribution of artificial SNe la is first 
resampled in z, s, c, and Amag. The goal is to represent the underlying population of (non- 
peculiar) SNe before selection effects by adopting distributions of light-curve width, color, 
and (Tint that match the observed distributions as a function of redshift in the selected sample. 
The low-redshift portion of the sample {z < 0.6) is extremely useful for constraining these 
distributions, but note that we are making the assumption that the underlying demographics 
of the population do not change significantly more strongly than the A + B model of Eq. [2] 
suggests. 

Rather than re-sample the input distributions, we instead simply weight the distribu- 
tions to obtain the same effect while obtaining higher statistics in our measurement by not 
throwing out an y objects. We wei g ht the SN distributio n to (Tint = 0.11 mag based on the cos- 



mological fits in lGuy et al.l ( l2010l ): IConley et al.l (120101 ). This differs from the value adopted 



therein because the simulated SNe used here do not include uncertainties in the input light- 
curve model or other systematics which tend to increase the measurement uncertainties. We 
find that the underlying stretch distribution is well described by a two-component Gaussian 
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distribution. A two-component Gaussian could not provide a good match to the observed 
color distribution and the evolution in the observed mean color with redshift due to selec- 
tion effects. Instead, we have adop ted the convoluti on of a exponential distribution and a 
Gaussian, similar in form to that of iJha et al.l (120071 ). 

The parameters of the stretch model are (mean, cr, peak value) 0.92,0.11,0.37 for the 
first Gaussian and 1.06, 0.08, 0.80 for the second. For the color model, the o of the Gaussian is 
0.04, the r of the exponential distribution 0.155, and the mean value of the distribution 0.07. 
Note that we do not interpret the color distribution as arising from host-galaxy extinction, 
but simply use this model as a useful empirical parameterization. The resulting color and 
stretch distributions are compared with the actual 3rd-year cosmological sample in Figure [TJ 
and the evolution of the mean parameters with redshift are shown in Figure El 



4.2. Recovery results 

The SN la recovery results are presented as a function of various parameters in Figures[9]- 
[TTl These plots include the cumulative data for all realizations performed on the field-seasons 
up to the end of D3 observations in June 2007. In each of the top panels of the figures, 
the input (resampled) distributions are shown as the unfilled histograms and the recovered 
populations are shaded in the lightest grey. The recovery data naturally include the timing, 
spacing, and quality of the actual observations, as we had the opportunity to detect each 
supernova while it was visible on any night that SNLS obtained usable images. 

The medium-grey histograms show the effects of incorporating a set of observational 
constraints on the artificial SNe that are recovered. These cuts ar e for the purpose o f 



ensuring adequate light-curve coverage for fitting, as also described in lPerrett et al.l (120 lOl ). 
The observational requirements are established as: 

• > 1 early (— 15d to +2.5d) observation in each of iu and r^; 

• > 1 early (— 15d to +5d) observation in gu] 

• > 1 observation near peak (— 9d to +7d) in each of iu and r^; 

• > 1 late (+5d to +20d) observation in either %m or tm- 

These dates are in effective rest-frame days, corrected for stretch and redshift: 

agGobs / . N 
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Similar cuts have been made in determining reliable efficiencies in other SN surveys (e.g. 



Dildav et al.ll2008l . l2010f ). 



The lowest, dark ffiled histograms in each top panel of Figures [UHTT] all represent the 
spectroscopically-viable population. In this case, the peak brightness and %inc limits are 
applied as described previously in ^ for follow-up suitability. In addition, the candidate 
must have been discovered prior to +7 days past peak (rest frame), and have at least one 
observation above background both before and after maximum light. Due to the limited 
availability of telescope time, not all of the suitable candidates would have actually been 
observed spectroscopically. 

Figure |9] shows the recovery fraction as a function of iM peak magnitude. The ability 
of the SNLS real-time pipeline to recover SNe la begins to drop off at magnitudes fainter 
than iM ~ 23 mag, with a 50% recovery limit at zab = 24.3 mag. This is measured in one 
individual night of observation; deeper limits can be achieved by combining multiple epochs 
of data in offline searches at the end of a run. The dip in the recovery fraction seen after 
the application of the observational cuts is likely due to the longer period over which the 
higher-stretch bright SNe la can be observed. 

The redshift distribution for the recovered objects (Figure ITOl-left) is fiat out to 2 ~ 0.5 
and falls off smoothly to a 50% incompleteness at ^ ~ 1. The likelihood of recovering SNe la 
is very small at 2; > 1.1 using this detection method. The small rise at low- 2; in the recovered 
fraction once the observational cuts are applied is due to time-dilation effects on the SN light 
curves (Eq. H]). Applying the observational constraints and spectroscopic cuts brings the 50% 
limit down to z ~ 0.85, as shown by the dotted line in the lower- left panel of Figure [TOl 

These results appear to b e consistent with si milar simulations carried out for the lower- 



redshift SDSS-II SN Survey (IDildav et all 120081 ). The SDSS search to z ~ 0.4 with an 



efflciency of 0.5 (compared with z ~ 1 for the SNLS, a factor of 2.5 larger), and report 
essentially 100% efflciency to 2 ~ 0.2 (compared with z ~ 0.5 in our simulations). The 
efflciencies in the SDSS are in good agreement with those from the SNLS, but multiplied by 
a factor 0.4 on the redshift axis. 

The recovery fraction as a function of stretch (Figure ITOl-center) exhibits an overall rise 
for s > 0.8 which is primarily a color effect: on average, lower-stretch SNe are redder, fainter, 
and tend to have smaller recovery efflciencies. The dip at s = 0.8 is an artifact of the break in 
the stretch-color relationship of Eq. |31 As expected, the recovery rate shows a clear decrease 
towards redder colors (Figure ITUI-right). 

The efflciency of SN la discovery depends not only on the brightness of the candidates, 
but also on the quality of the background subtraction in the real-time images. Furthermore, 
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candidates with low %inc values (Eq. [T]) compared with their hosts will generally not be con- 
sidered suitable for spectroscopy. Figure [TT] provides the distributions and recovery fractions 
as a function of Amag (left), galactocentric radius (center), and %inc (right). The Amag 
recovery shows a modest decrease towards increasing values due to the reduced chance of 
finding intrinsically fainter objects. The recovery fraction is constant for positions out to 
r ~ 1.0" from the host center, beyond which radius the relative effects of any galaxy sub- 
traction residuals will begin to diminish and the detection efficiency will increase. Note that 
these effects are averaged over hosts of various morphological types and brightnesses; the 
PSF convolution and subtraction of more highly structured galaxies can leave residuals in the 
difference images that can conceal detections. The %inc plots at the right of Figure [TT] reveal 
that the fractional brightness over the host galaxy background can have a significant impact 
on detection recovery, although the effects of candidate brightness and color are also folded 
into this effect. The details of the host selection process in the artificial SN la simulations 
(Eq. [2]) are found to have little effect on the recovery results. 

Some of the recovery plots in Figures [TOl and [TT] can be difficult to interpret in view of the 
fact that they incorporate a broad and uniform distribution of redshifts from z = 0.1 — 1.2. 
To address this. Figure ITU provides the spectroscopic recovery fractions in four redshift bins 
as a function of Amag, c, host radius, and %inc. In each case, the incompleteness in the 
high-redshift bins is the most severe; this is due primarily to the brightness and %inc limits 
imposed on the spectroscopic sample. 

The human review process is another element that must be considered in the detection 
recovery estimates; this is complicated by the fact that fake SNe were not sho wn to human 



reviewers during the RTA, as was done in some surveys (e.g.. lDilday et al.ll2008l ). In the RTA 
procedure, candidates were visually confirmed by an operator who could remove obviously 
erroneous detections prior to updating the SNLS database. Since any false detections that 
made it past the reviewer into the database could later be culled based on measurements 
of additional photometry and subsequent light-curve fits, the only remaining concern is the 
omission of good SN la candidates due to human error. 



Neill et al.l (120061 ) performed a Monte Carlo simulation to test the SNLS human review 
process using single-epoch images having image quality (IQ) values of 0'.'69 and 1'.'09. Their 
results indicate that the Canadian RTA procedure yields a raw recovery percentage better 
than ~ 95% brightward of the incompleteness limit. Given multiple chances of candidate 
discovery using several epochs of SNLS data, this success rate has been found to approach 
~ 100% over the duration of a typical, well-observed SNe la. Therefore, any additional 
contribution of the visual review process to candidate recovery is deemed to be insignificant. 

A comparison of the input and measured photometry for the artificial SNe la in the 



- 26 - 



simulations is shown in Figure [131 No significant systema tic biases a r e foun d in the measured 
photometry using the RTA pipeline. In bright galaxies, iGuy et al.l (120101 ) have identified a 
~ 2% oversubtraction of host galaxy flux at the location of the SN when applying a similar 
method in the flnal photometry analysis; this bias is not apparent to within the precision of 
the RTA measurements. 



5. Selection bias in the SNLS sample 



The effects of systematic errors introduced by Malmquist bias and spectroscopic sam- 
ple selection ( §3.ip are important considerations for supernova cosmology: they cause a 
redshift-dependent offset in the average relative distance moduli from which the cosmologi- 
cal parameters are determined. The distance moduli (/i^) for a Type la SN can be calculated 
from: 

i|j - M + a(s - 1) - /3c, (5) 



m. 



where m*^ is the peak magnitude in the B band, s and c are stretch and color. The absolute 
magnitude, M, and t he coefficients a an d /3 are determined by minimizing the residuals on 
the Hubble Diagram (lAstier et al.ll2006l ). 



The artificial SNe la from §4.21 are used to calculate the mean offsets in Amag, stretch, 
and color as a function of increasing redshift. The mean properties of the observed spectro- 
scopic sample can be compared with the underlying (unbiased) population to look for the 
effects of sampling biases. The unbiased population is represented by the input distribu- 
tions that were re-weighted to match the underlying distribution as discussed in §4.11 (open 
histograms in the upper panels of Figures lU HTT]) . This, of course, assumes that there is no 
significant evolution in the intrinsic properties of SNe la with redshift. 

The observed spectroscopic sample consists of recovered objects that were deemed suit- 
able for follow-up. We also assume that any SNLS candidates that were observed spec- 
troscopically could be randomly selected from this latter sample; due to observing time 
constraints, not all suitable objects were queued for spectra. Furthermore, we do not include 
a component to address the spectroscopic identification efficiency. This would take into 
account the candidates that were observed spectroscopically, but which were unidentifiable 
or misclassified from their spectra. Such sample contamination issues are expected to be 
relatively min or for bias calculations, yet are exceediri gly difficult to determine precisely (see 
discussions in iBalland et al.l l2009l : IConley et al.l l2010l ) . 



The offset in the mean Amag for the observed sample as a function of redshift is pre- 
sented in Figure [TH showing the combined effects of Malmquist bias and spectroscopic 
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Table 4. The SNLS sample bias for the magnitude dispersion (Amag), SN stretch (s) and 

SN color (c). 



z 


(^(Ama 


g) 




















0, 


.33 


+0.001 ± 


0, 


.001 


0, 


.007 + 0, 


c\r\c\ 

.002 


— 0, 


r\c\-\ 

.001 


1 

+ 


0, 


c\r\ 1 

.001 


0. 


,37 


+0.000 ± 


0, 


.001 


0, 


.006 + 0, 


.001 


— 


.002 


+ 


0, 


.001 


0, 


,42 


—0.001 ± 





.001 


0, 


.006 + 0, 


.001 


— 


.002 


+ 


0, 


.000 


0. 


,47 


—0.002 ± 


0, 


.002 


0, 


.007 + 0, 


.002 


— 


.006 


+ 


0, 


.002 


0, 


,53 


-0.003 ± 





.002 


0, 


.009 + 0, 


.002 


-0 


.009 


+ 


0, 


.002 


0. 


,58 


-0.004 + 


0, 


.001 


0, 


.009 + 0, 


.003 


-0, 


.014 


+ 


0, 


.003 


0. 


,62 


-0.007 + 


0, 


.001 


0, 


.011 + 0, 


.002 


-0, 


.019 


+ 


0, 


.002 


0. 


,68 


-0.009 + 


0, 


.002 


0, 


.015 + 0, 


.003 


-0, 


.027 


+ 


0, 


.003 


0. 


,73 
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.001 
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.004 
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+ 
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.053 
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0, 
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+ 


0, 
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0. 
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-0.024 + 





.001 


0, 
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.004 


-0 


.073 


+ 


0, 


.003 


1. 
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-0.029 + 


0, 


.004 


0, 
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.002 
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+ 


0, 


.002 


1. 
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-0.038 + 


0, 


.003 


0, 


.053 + 0, 


.005 
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.002 
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,13 
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0, 


.003 
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.110 
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^For the parameters Amag, s and c, the numbers listed 
in each redshift bin arc the difference, (5, between the 
mean of each parameter for the input SN population, and 
the mean of the same parameter for the recovered SN 
sample. The errors quoted are the errors in the mean. 
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sampling. The sampling effects start to become significant at z > 0.75, but are apparent 
starting at 2 ~ 0.6, shifting the mean Amag in the observed sample towards brighter values. 
Column 2 of Table H] provides the numerical values of the average offsets in redshift bins for 
all fields combined. The mean magnitude corrections measured for the individual fields are 
all consistent and are listed in Table [51 

The offset in observed Amag with redshift depends strongly on the intrinsic dispersion of 
the actual sample. A comparison of the mean bias for different assumed values of a"int in the 
artificial population is given in Figure [T^ The lower the dispersion in the Amag distribution, 
the smaller the effects of selection bias. This is due to a reduction in the relative contribution 
from intrinsically faint input objects near the detection limit. 

The stretch and color trends with redshift were shown in Figure [Hi The increasing offset 
towards larger stretch results from preferentially finding the brighter objects with broader 
light curves at high z (Column 3 in Table [1]). Similarly, the mean color offsets (Column 4 
in Table [1]) are caused by a bias towards detecting brighter — and hence bluer — objects 
at 2; > 0.8. These stretch and color offsets are not directly relevant to the cosmology: we 
already correct for first-order variations in s and c via the a and /? parameters in EqJ^l and 
the second-order sampling biases that affect the average intrinsic luminosity of SNe la with 
redshift are already incorporated into the Amag corrections provided above. 



6. Systematic Errors in the Malmquist Bias 

There are many parameters in our estimate of the Malmquist bias, including: the as- 
sumed color-color law, the parameters of the spectroscopic selection model and its form, the 
values of cxint, a, /3 and M, and the measured underlying stretch and color distributions. In 
this section, we discuss the effects of these terms on the Malmquist bias. 

As noted previously, the color-color relation used to warp our SN template to match 
the simulated colors does not match our current understanding of SN behavior. In order to 
determine the effects of this on our r esults, we re-weig hted the simulations to reproduce the 



effects of using the SALT2 color law (iGuy et al.ll2007l ). The effects are quite small, less than 



0.0003 mag at 2; = 1, so we have not applied this correction to the results. 

The uncertainties in the luminosity distance model are 0.01 mag for cjint, 0.08 for a, 0.1 
for /3, and 0.025 for Ai0. We assign uncertainties of 0.1 mag to each of the magnitude limits 



^^We have assumed Hq = 70 km/sec/Mpc. Formally, our distance model is only sensitive to some combi- 
nation of Ho and M, and it is this quantity that is constrained to 0.025 mag, independent of any assumptions 
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Table 5. The SNLS sample bias on the magnitude dispersion (Amag) in each search field. 



z 


5(Amag)Di'' 


(^(Amag)D2 


5(Amag)D3 


5(Amag)D4 


U.oo 


1 n nno 
+U.UUZ 


1 n nni 
+U.UUi 


—U.UUi 


1 n nm 
+U.UUi 


U.d t 


— U.UUi 


+U.UUi 


1 n nm 
+U.UUi 


1 n nnn 
+U.UUU 


n /I o 
U.4z 


—U.UUi 


n nm 
—U.UUi 


n nno 
— U.UUz 


n nnn 
—U.UUU 




— U.UUU 


n nn/1 
— U.UU4 


n nn i 
—U.UUi 


n nn i 
—U.UUi 


0.53 


-0.001 


-0.005 


-0.003 


-0.004 


0.58 


-0.004 


-0.004 


-0.003 


-0.006 


0.62 


-0.006 


-0.009 


-0.006 


-0.006 


0.68 


-0.008 


-0.012 


-0.007 


-0.010 


0.73 


-0.009 


-0.007 


-0.009 


-0.010 


0.77 


-0.007 


-0.011 


-0.012 


-0.010 


0.82 


-0.013 


-0.012 


-0.012 


-0.015 


0.88 


-0.012 


-0.018 


-0.015 


-0.016 


0.92 


-0.017 


-0.022 


-0.016 


-0.015 


0.97 


-0.024 


-0.024 


-0.023 


-0.023 


1.03 


-0.023 


-0.032 


-0.032 


-0.030 


1.07 


-0.038 


-0.041 


-0.034 


-0.040 


1.13 


-0.052 


-0.053 


-0.054 


-0.059 


1.17 


-0.072 


-0.084 


-0.066 


-0.059 



^In each redshift bin, the difference, 5, between the mean 
Amag for the input SN popiilation and the mean of the same pa- 
rameter for the recovered SN sample is given. The errors quoted 
are the errors in the mean. 
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of our spectroscopic selection model, and 25% to each of the %inc components. In addition, 
we consider changing the form of the spectroscopic model to linearly interpolate between 
the %inc limits as a function of magnitude rather than abruptly changing from one limit to 
another. For the color and stretch models, we vary the parameters of the distributions until 
the of the fit to the stretch and color evolution with redshift increases by one. The effects 
of each of these components are given i n Table [6] in term s of the mean Malmquist bias of 
the 3rd- year SNLS cosmo logical sample ( iGuy et al.l 120101 ) and the bias at 2; = 1. Since the 
mean and z = 1 values are 0.0065 and 0.0237 mag, respectively, the total systematic error is 
approximately 20%. 



7. Summary 



The SNLS real-time analysis has discovered thousands of supernova candidates during 
its first four years of full operation. The vast majority of these candidates remain uncon- 
firmed by spectroscopy due to observing time limitations, although methods for photometric 
identification are improving. Nonetheless, this will remain a major limitation for the next 
generation of surveys that could discover an order of magni tude more SNe, yet may not have 
access to much more in the way of spectroscopic time (see lHowell et al.l 120091 ). 



Supernova surveys all suffer from the effects of selection biases that must be carefully 
considered in any analysis (e.g., iKessler et al.ll2009l ). Monte Carlo simulations of SNe la 
in the SNLS real-time pipeline show that the effects of Malmquist bias and spectroscopic 
sample selection begin to affect the distribution of intrinsic magnitudes observed at 2; > 0.6. 
At 2; = 1, the offset in the mean Amag reaches ~ —0.027 mag as the survey systematically 
observes the intrinsically brighter SNe. We find approximately 20% systematic errors in 
our estimate of the Malmquist bias based on the uncertainty in the input model to our 
simulations. 



z ~ 



The SNLS has produced a supernova dataset of unprecedented size and quality out t o 
1.1. This large, homoge neous sample of SNe la is ideal for cosmoloev (lAstier et al. I l2006h 



the analysis of SN la rates (iNeill et al. 



2006 



Perrett et al. 



2008 



2010 



Carlberg et al. 



20081: 



Ripoche et aPboioh. and 



Sullivan et al. 



studi es of the SN environments ( iGraham et al. 
20101). Other types of investigations can be made usi ng the various sarn ples of objects 
in the SN LS database: e.g ., cosmology with SNe IIp (Nugent et al. 20061), co re-collapse 
SN rates ( Bazin et al. 2009), sub-lu r ainous SNe la JConzalez-Gaitan et al. I boioh . analyses 



of unusual SNe ( iHowell et al.l 120061 : ISchawinski et al.ll2008l ). as well as optical variability 



about Hq. 
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Fig. 7. — The recovered stretch (left) and color (rig ht) distributions compared with those 
observed for the 3rd-year SNLS cosmo logical sample ( iGuy et al.ll2010l ). 



Table 6. Systematic uncertainties in the Malmquist bias estimate arising from 

assumptions in our simulations. 



Term 


Average systematic 


Total uncertainty 




uncertainty 


at 2; = 1 


Spec model params 


0.0006 


0.0015 


Spec model form 


0.0002 


0.0010 


C"int 


0.0013 


0.0049 


a 


0.0002 


0.0005 




0.0003 


0.0008 


M 


0.0002 


0.0003 


Stretch Distribution 


< 0.0001 


0.0002 


Color Distribution 


0.0002 


0.0004 


Color-color law 


< 0.0001 


0.0003 


Total 


0.0015 


0.0053 
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Fig. 8. — The recovered stretch evolution in the mean stretch and color (blue squares) 
compared with the actual stretch and color distributions (averages shown as red circles) as 
a function of redshift. 
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Fig. 9. — The upper panel shows the input (unfilled) and recovered (light grey) distributions 
in peak zm magnitude for the artificial SNe la used in the recovery tests. The efi^ects of adding 
the observational and spectroscopic constraints are shown by medium and darkest grey 
histograms, respectively. The lower panels present the recovered fraction for the found sample 
(solid line) and with the observational and spectroscopic cuts (dashed and dotted lines, 
respectively). The 50% incompleteness limit for the real-time analysis lies at zab = 24.3 mag 
(AB), represented by the vertical line in the lower diagram. 
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Fig. 10. — SN la recovery as a function of redshift (left), stretch (center), and color (right). 
The histograms in the upper panels show the input distributions (unfilled) and the recovered 
sample (light grey), along with the addition of the observational and spectroscopic cuts 
(medium and darkest grey histograms, respectively). The lower panels present the recovered 
fraction for the found sample (solid line) and with the observational and spectroscopic cuts 
(dashed and dotted fines, respectively). 
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Fig. 11. — SN la recovery as a function of Amag (left), galactocentric radius in the host (cen- 
ter), and %inc (right), as in Figure [TUl Note that the detection effects are averaged over hosts 
with a variety of morphological types. The bar at the left of the %inc histogram represents 
all objects with %inc< 0, while the bar at the right includes SNe la with Log(%inc)> 4. The 
large variations in the recovery fraction at Log(%inc)< 1 are due to small sampling statistics 
in those bins. 



- 36 - 



100 



100 




0.5 1.0 1.5 
Host Radius (arcsec) 



1 2 3 

Log(% Increase) 



Fig. 12. — SN la recovery fraction for the spectroscopic sample shown separately in four 
redshift bins, given as a function of Amag, color, galactocentric radius in the host, and 
%inc. The lines represent the following redshift ranges: 0.10 < z < 0.35 (solid blue), 
0.35 < z < 0.60 (long-dashed green), 0.60 < z < 0.85 (short-dashed orange), 0.85 < z < 1.10 
(dotted red). In each case, the high-redshift range (dotted red line) has the lowest recovery 
rate, due mostly to the peak magnitude and %inc restrictions on the spectroscopic sample. 
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Fig. 13. — Contours showing the difference between input and measured magnitude (zm at 
peak) versus input magnitude for the recovered sample of artificial SN la (top), and those 
missed by the RTA detection pipeline (bottom). The vertical lines show the 50% detection 
incompleteness hmit of im — 24.3 mag (AB). 
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Fig. 14. — Malmquist and spectroscopic selection bias in the SNLS sample, as shown by the 
mean offset in Amag as a function of redshift. The errors represent simply the statistical 
uncertainty from the simulations. The small error aX z = 0.97 is a statistical fluctuation. 
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Fig. 15. — The mean offset in Amag as a function of redsliift, as in Figure [HI shown for 
three different values of cxinf The actual uncertainty in ai-at is around 0.01 mag. 
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studies of AGNs. In a carefully designed survey, the process of finding supernovae is quite 
straightforward. The real difficulty lies in determining accurate calibrations and in the 



management of detailed s ystematic errors (lAstier et aLKOOGuRegnault et al.ll2009l : iGuy et al. 



2OIOI : IConley et al.ll2010l ). Future supernova surveys at all redshifts must focus on obtaining 



more high-quality and well-characterized data, not simply on detecting more supernovae. 
That challenge will demand significantly more in the way of telescope time and resources. 
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